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Abstract 

A striking aspect of the recently proposed split supersymmetry is the existence of heavy gluinos 
which are metastable because of the very heavy squarks which mediate their decay. In this paper 
we correlate the expected flux of these particles with the accompanying neutrino flux produced in 
inelastic pp collisions in distant astrophysical sources. We show that an event rate at the Pierre 
Auger Observatory of approximately 1 yr^^ for gluino masses of about 500 GeV is consistent with 
existing limits on neutrino fluxes. Such an event rate requires powerful cosmic ray engines able 
to accelerate particles up to extreme energies, somewhat above 5 x 10^^ GeV. The extremely low 
inelasticity of the gluino-containing hadrons in their collisions with the air molecules makes possible 
a distinct characterization of the showers induced in the atmosphere. Should such anomalous events 
be observed, we show that their cosmogenic origin, in concert with the requirement that they reach 
the Earth before decay, leads to a lower bound on their proper lifetime of the order of 100 years, and 
consequently, to a lower bound on the scale of supersymmetry breaking, AsusY > 2.6 x 10^^ GeV. 
Obtaining such a bound is not possible in collider experiments. 

PACS: ll.30.Pb, 96.40.-z, 13.85.Tp 
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I. GENERAL IDEA 



The standard model (SM) of particle physics has had outstanding success in describing all 
physical phenomena up to energies ~ 500 GeV [1]. Nonetheless, there is a general consensus 
on that it is not a fundamental theory of nature (apart from the fact that the SM docs not 
include gravity): with no new physics between the energy scale of clectroweak unification 
{Mw ~ 10^ GeV) and the vicinity of the Planck mass (Mpi ~ 10^^ GeV) the higgs mass must 
be fine-tuned to an accuracy of order (Mvy/Mpi)^ to accommodate this enormous desert. The 
leading contender for the elaboration of the desert has been the supersymmetric extension 
of the SM [2]. Supersymmetry (SUSY) posits a "complete democracy" between integral and 
half-integral spins, implying the existence of many as- yet-undiscovered superpartners. Thus, 
if SUSY can serve as a theory of low energy interactions, it must be a broken symmetry. The 
most common assumption is that the minimal low energy effective supersymmetric theory 
(MSSM) has a breaking scale of order Asusy ~ 1 TeV, thus avoiding 't Hooft naturalness 
problem with the higgs mass. 

MSSM has a concrete advantage in embedding the SM in a grand unified theory: the 
supersymmetric beta functions for extrapolating the measured strengths of the strong, elec- 
tromagnetic, and weak couplings lead to convergence at a unified energy value of the order 
Mgut ~ 10^^ GeV [3]. The model, however, is not free of problems. In particular, dimension 
four i?-parity violating couplings in the superpotential yield unacceptably large proton decay 
rates and neutrino masses. This can be readily solved by imposing i?-parity conservation, 
which as a byproduct ensures the stabihty of the lightest SUSY particle, making it a possible 
candidate for cold dark matter [4]. However, there arc other problems in the MSSM: di- 
mension five operators tending to generate excessive proton decay; new CP-violating phases 
which require suppression for agreement with limits on electric dipole moment limits; and 
excessive flavor violations, due e.g., to the absence of a complete flavor degeneracy in the 
Kahler potential of minimal supergravity [5] . 

Of course, the fine-tuning involved in accommodating the above constraints is miniscule 
in comparison to that required in generating a cosmological constant that satisfies 't Hooft 
naturalness. Recent experimental data [6] strongly indicate that the universe is expanding in 
an accelerating phase, with an effective de Sitter constant H that nearly saturates the upper 
bound given by the present-day value of the Hubble constant, i.e., H ^ Hq ^ 10~^^ eV. 
According to the Einstein field equations, provides a measure of the scalar curvature 
of the space and is related to the vacuum energy density, £4, according to M|ji7^ ~ £4. 
However, the "natural" value of £4 coming from the zero-point energies of known elementary 
particles is found to be at least €4 ~ AgusY; yielding H > 10~^ eV. The failure of 't Hooft 
naturalness then centers on the following question: why is the vacuum energy determined by 
the Einstein field equations 60 orders of magnitude smaller than any "natural" cut-off scale 
in effective field theory of particle interactions, but not zero? Nowadays, the only existing 
framework which can address aspects of this question is the anthropic approach [7] . 

Very recently Arkani-Hamed and Dimopoulos (AD) have looked at SUSY from a different 
angle [8]. In their model, the scale of SUSY breaking is pushed to a very high energy (say, 
Asusy ~ 10^"^ GeV) and a higgs mass of order TeV is recovered by invoking fine tuning. For 
this breaking scale, the bosonic superpartners are heavy, while the extra fermions retain TeV- 
scale masses thanks to protection by chiral symmetry. (We follow Giudice and Romanino [9] 
in adopting the designation "split SUSY" for the AD model). This scenario preserves the 
achievements of the MSSM while resolving the problems mentioned above. In particular. 
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analyses of one loop [10] and two loops [9] running of the RG equations, show that the 
AD scenario preserves unification of couplings. Moreover, aside from the light higgs tuning, 
the other flaws inherent to the MSSM elegantly disappear when the scalar superpartners 
decouple. 

The AD model can be discussed in the same anthropic framework adopted for exami- 
nation of the cosmological constant. Recent investigations in String Theory have applied 
a statistical approach to the enormous "landscape" of vacua present in the theory [11]. 
Among this vast number of metastable vacua, there can be small subset C(10^°) exhibiting 
low scale SUSY breaking, a TeV-scale higgs, as well as the remaining traditional MSSM 
physics [12]. However, the fine tuning required to achieve a small cosmological constant im- 
plies the need of a huge number of vacua, far more than the (9(10^°) characterizing low-scale 
SUSY breaking [13]. Remarkably, if one posits high-scale SUSY breaking and superpartners 
widely separated in mass, 0(10'^^^) vacua become available, enough to fine tune both the 
cosmological constant and the higgs mass. 

It is therefore instructive to explore how drastically the AD scenario can change the 
phenomenology of conventional MSSM. Prospects for probing split SUSY at the LHC [14] 
as well as in dark matter searches [15] have been recently developed. In what follows we 
show that cosmic ray data may also provide important information about the AD scenario. 
As a principal result of this paper, we will delineate conditions under which one can set a 
lower bound on the SUSY breaking scale in a region of parameter space far beyond that 
probed at the LHC. 

An intriguing prediction in this scenario, which represents a radical departure from the 
MSSM, is the longevity of the gluino [16]. As mentioned above, in split SUSY the squarks are 
very massive and so gluino decay via virtual squarks becomes strongly suppressed, yielding 
a g lifetime of the order of [8] 



where Mg is the gluino mass. Very strong limits on heavy isotope abundance in turn re- 
quire the gluino to decay on Gyr time scales, leading to an upper bound for the scale of 
SUSY breaking 0(10^^) GeV [8]. Because of the large mass of the gluino, the threshold for 
inelastic scattering on the cosmic microwave background is ^ 10^"^ GeV [17], allowing ultra- 
high energy gluino-containing hadrons ("G's") to reach us unimpeded from cosmological 
distances. 

In this work we study the possibility of G-detection with cosmic ray observatories. To 
this end, in Sec. II we discuss the main characteristics of cascades induced by G-hadrons and 
estimate the sensitivity of the Pierre Auger Observatory (PAO). After that, in Sec. HI, we 
correlate the expected fiux of G-hadrons with the accompanying neutrino fiux produced in 
inelastic pp collisions in distant astrophysical sources. We show that an event rate ~ 1 yr~^ 
at PAO for gluino masses of about 500 GeV is consistent with existing limits on neutrino 
fluxes. The actual observation of a few G-cvents will then directly imply a lower bound 
on AsusY- The details of this interesting possibility are presented in Sec. IV. Section V 
contains our conclusions. 




(1) 



3 



II. CHARACTERISTICS OF AIR SHOWERS INITIATED BY G-HADRONS 



The interaction of a high energy cosmic ray in the upper atmosphere gives rise to a roughly 
conical cascade of particles that reaches the Earth's surface in the form of a giant "saucer" , 
traveling at nearly the speed of light. In the case of proton- or nucleus-induced cascades, 
the leading particle and other high energy hadrons (mostly pions) in the shower core readily 
cascade to lower energies as they interact with the air molecules. Because of the prompt 
decay of neutral pions, 1/3 of the energy in each of these interactions transits into energetic 
7-rays. Electromagnetic subshowers are then initiated: the high energy photons produce 
pairs that lose energy by bremsstrahlung and ionization before annihilation into a new 
photon at lower energy. Eventually, the average energy per particle drops below a critical 
energy eo ~ 86 MeV at which point ionization takes over from bremsstrahlung and pair 
production as the dominant energy loss mechanism. The changeover from radiation losses 
to ionization losses depopulates the shower, and defines Xmax, the longitudinal coordinate 
of maximum multiplicity. 

The number of muons (and neutrinos) does not increase linearly with energy, because at 
higher energy more generations are required to cool the pions to the point where they are 
likely to decay before interaction. Production of extra generations results in a larger fraction 
of the energy being lost to the electromagnetic cascade, and hence a smaller fraction of the 
original energy being delivered to the tt^. The electrons, positrons and photons are thus the 
most prolific constituents in the thin disk of particles showering towards the ground, and 
most of the energy (about 90%) is dissipated in the electromagnetic cascade. 

By the time they reach the ground, relatively vertical showers have evolved fronts with 
a curvature radius of a few km, and far from the shower core their constituent particles are 
well spread over time, typically of the order of a few microseconds. For such a shower both 
the muon component and a large portion of the electromagnetic component survive to reach 
the ground. For inchned showers the electromagnetic component is absorbed long before 
reaching the ground, as it has passed through the equivalent of several vertical atmospheres: 
2 at a zenith angle 9 = 60°, 3 at 70°, and 6 at 80°. In these showers, only high energy muons 
created in the first few generations of particles survive past 2 equivalent vertical atmospheres. 
The rate of energy attenuation for muons is much smaller than it is for electrons, thus the 
shape of the resulting shower front is very fiat (with curvature radius above 100 km), and 
its time extension is very short (less than 50 ns) [18]. 

The energy lost by a G during collision with nucleons is primarily through hard scatter- 
ing [19]. This imphes a fractional energy loss per coUision, i^Tinei ~ 1 GeV/Mo- A heuristic 
justification for this result is as follows: consider the inclusive process GN GX. It is a 
simple kinematic exercise to show that the minimum momentum transfer for Mx, Mq <^ s 
is given by \t,ni,i\ — Mq Mx/s^^ where ^/s is the center-of-mass energy of the GN collision. 
Kinematics also determine the fractional energy loss, i^Tmei = M\j s. Combining these two 
results we obtain that Xinei — |^mmP''^/-^G- The hard scattering restriction in QCD requires 
\imin\ ~ 1 GeV^, so that for a fiducial Mq = 500 GeV, we obtain Xinei ~ 0.002. Note that, 
for Mg = 50 GeV, our formula agrees with the result in [19]. 

In the case of a G-induced shower, the very low inelasticity of G-air interactions implies 
the leading particle retains most of its energy all the way to the ground, while the secondary 
particles promptly cascade to low energies as for any other air shower. This results in 
an ensemble of mini-showers strung along the trajectory of the leading particle. Since 
the typical distance between mini-showers is about 10 times smaller than the extent of a 
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single longitudinal profile, it is not possible to resolve the individual showers experimentally. 
Instead one observes a smooth envelope encompassing all the mini-showers, which extends 
from the first interaction all the way to the ground. Monte Carlo simulations have been 
performed which confirm this phenomenological description (see Fig. 3 in Ref. [19]). The 
G-showers indeed present a distinct profile: (1) there is only a few percent probability for 
its Xjnax to be mistaken for that of a proton shower [20] (2) the flatness of the longitudinal 
development is unique to the extremely low inelasticity of the scattering, and can be easily 
isolated from background. 

The characteristics of G-induced showers observed at the ground should also be distinct 
from those characteristic of proton and nucleus induced showers. Each mini-shower gener- 
ates a bundle of muons which survive to the ground. Since each bundle is produced at a 
different slant-depth, the muon component of the shower front exhibits a much more pro- 
nounced curvature than what one would expect for a shower of the same energy interacting 
only in the upper atmosphere. The difference between a proton/nucleus and a G would thus 
be much more evident in the case of inclined showers for which there is much less electromag- 
netic contamination. Specifically, a highly inclined G-induced shower produces many muon 
bundles, and so should exhibit a much more curved shower front than a proton/nucleus 
shower with the same energy and zenith angle. 

The experimentally interesting region to search for G events then lies between 70° < 
9 < 90°. The reduction of the solid angle acceptance to larger than 70° eliminates the 
hadronic background. Moreover, there needs to be sufficient pathlength for the G, with 
its low inelasticity, to lose sufficient energy. The mean free path for ultra-high energy G's 
is about 30 g cm~^ [19], and the pathlength of the atmosphere at 70° is about 100 times 
this [18], allowing more than 20% of the energy to evolve in the air shower. 

With this in mind, the G-event rate for a given cosmic ray experiment is found to be 



where A{Eg) is the hadronic aperture for 9 > 70°. 

There are two major techniques which can be employed in detecting cosmic ray air show- 
ers. The most commonly used detection method involves sampling the shower front at a 
given altitude using an array of sensors spread over a large area. Sensors, such as plastic 
scintillators or Cerenkov detectors are used to infer the particle density and the relative 
arrival times of the shower front at different locations. The muon content is usually sought 
either by exploiting the signal timing in the surface sensors or by employing dedicated detec- 
tors which are shielded from the electromagnetic shower component. Inferring the primary 
energy from energy deposits at the ground is not completely straightforward, and involves 
proper modeling of both the detector response and the physics of the first few cascade gen- 
erations. Another highly successful air shower detection method involves measurement of 
the longitudinal development of the cascade by sensing the fluorescence light produced via 
interactions of the charged particles in the atmosphere. Excited nitrogen molecules fluo- 
resce producing radiation in the 300 - 400 nm ultraviolet range, to which the atmosphere 
is quite transparent. The shower development appears as a rapidly moving spot of light 
whose angular motion depends on both the distance and the orientation of the shower axis. 
The integrated light signal is proportional to the total energy deposited in the atmosphere. 
Fluorescence observations can only be made on clear moonless nights, yielding a duty cycle 
of about 10%. 




(2) 
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Over the next few years, the best observations of the extreme end of the cosmic ray spec- 
trum will be made by the PAO [21], which is actually comprised of two sub-observatories. 
The Southern site is currently operational and in the process of growing to its final size of 
S ~ 3000 km^. Another site is planned for the Northern hemisphere. The PAO works in a 
hybrid mode, and when complete, each site will contain 24 fluorescence detectors overlook- 
ing a ground array of 1600 water Cherenkov detectors. During clear, dark nights, events 
are simultaneously observed by fluorescence light and particle detectors, allowing powerful 
reconstruction and cross-calibration techniques. Simultaneous observations of showers us- 
ing two distinct detector methods will also help to control the systematic errors that have 
plagued cosmic ray experiments to date. 

For showers at inclination 6 > 70° and energy i?sh > 10^ GeV the probability of detecting 
an event falling within the physical area is roughly 1. Note that, from the previous discussion. 
Eg ^ 5 E'sh ^ 5 X 10^ GeV. The total aperture (2 sites) of the surface array is found to be 

r90° 

ApAo = S COS 9 dfl ^ 2200 km^ sr , (3) 

J70° 

where S cos 6 is the projected surface of the array in the shower plane, and dQ is the 
acceptance solid angle. From Eqs. (2) and (3), one notes that an event rate of ~ 1 yr~^ 
requires an integrated flux 

/ Jg(^g) dEo ^ 1.4 X 10-^^ cm-2 s"^ sr"^ . (4) 

J Eg, min 

Requiring that events also trigger the fluorescence detectors increases this flux by a factor 
of 10. 

Detailed characteristics of ultra-high energy cosmic ray sources are largely unknown [22] , 
rendering a direct calculation of the expected G-fiux speculative. In the next section we 
assess the viability of this flux by comparing with existing limits on gamma ray and neutrino 
fluxes. 



III. PRODUCTION OF G-HADRONS IN ASTROPHYSICAL SOURCES 

Among the non-thermal sources in the universe, radio-loud active galactic nuclei (AGNs) 
seem to be the most important energetically. There are other interesting powerful sources, 
like gamma ray bursts and cluster of galaxies; however, the non-thermal energy release in 
these astrophysical processes does not come close that of AGNs. At radio frequencies, where 
very large baseline interferometers can resolve the emission regions at milliarcsecond scale, 
many of radio-loud AGNs exhibit compact jets of relativistic plasma which are remarkably 
well coUimated, with opening angles about a few degrees or less. The AGNs come in various 
disguises according to the orientation of their radio jets axes and characteristics of the 
circum-nuclear matter in their host galaxies. The most extreme versions are Fanaroff Riley 
radio-galaxies with the radio jet axes almost in the plane of the sky and blazars with the 
radio jet axes pointing close to the line of sight to the observer, yielding a significant flux 
enhancement because of Doppler boosting. 

A total of 66 blazars have been detected to date as GeV 7-ray sources by the Energetic 
Gamma Ray Experiment Telescope (EGRET) on board of the Compton Gamma Ray Ob- 
servatory (CGRO) [23]. In addition, observations from ground-based Cerenkov telescopes 
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indicated that at least in 2 of these sources the 7-ray spectrum can be traced to more than 
a TeV observed photon energy [24]. The non-thermal emission of these powerful objects 
indicates a double-peak structure in the overall spectral energy distribution. The first com- 
ponent (from radio to X-rays) is generally interpreted as being due to synchrotron radiation 
from a population of non-thermal electrons, whereas the second component (7-rays) is ex- 
plained either through inverse Compton scattering of the same electron population with the 
various seed photon fields traversed by the jet [25] , or by the decay of neutral pions produced 
when the highly relativistic baryonic outflow coUides with diffuse gas targets moving across 
the jet [26]. 

In this work we focus on the "relativistic jet meets target" scenario, in which G-hadrons 
can be produced in collisions of ultra-high energy protons in the jet with those in surrounding 
gas. In the course of these collisions pions are produced which, on decay, give rise to a flux of 
photons and neutrinos. In what follows, we estimate the relative probabilities for production 
of G's and high energy neutrinos. We can then assess whether existing limits on the flux of 
high energy energy neutrinos and EGRET data from low energy gamma rays are consistent 
with the G flux given in Eq. (4). 

In order to specify detection criteria, the following kinematic analysis is relevant. The 
average energy of the produced G in the target system is given by 



Plab 

, (5) 



\ 2Mp 



where E^^^ is the energy of the high energy proton in the jet, and E(f^ is the G energy in 
the center-of-mass of the pp collision. Full acceptance at PAO requires a minimum energy 
-^c'min — 5 X 10^ GcV for the showering G hadron. From Eq. (5), this determines a minimum 
energy for the high energy proton: 



> 2 X 10» GeV (52i|!XV fiiffi 



where we have made use of the fact that, on the average, Eq^ = V^/2, where s is the square 
of the energy in the center-of-mass of the parton-parton collision producing the G's. 

It is immediately apparent from Eq. (6) that ultra-high energy sources are required in 
order to produce G's which can generate extensive air showers. The required energy can 
be reduced by restricting production to large s. In order that the maximum energy at the 
source does not exceed 10^^ GeV, a speculative number sometimes used in the literature [27], 
we take s > 16M^, which via Eq. (6) leads to E^^^^ 5 x 10^^ GeV. 

We now turn to evaluating the consistency of the required G flux and its accompanying 
pion flux with existing limits on neutrino and gamma ray fluxes. The required relationship 
is 

ME^) dE^ = I / JaiEa) dEa , (7) 



3 (Tpp^G{Smin) Nq 

where is the neutrino flux, cTinei and (Tpp^c are the total pp inelastic and inclusive pp ^ G 
cross sections, and N^, and Nq = 2 are the neutrino and gluino multiplicities per collision. 
The factor of 2/3 accounts for the fact that only charged pions contribute to neutrino 
production. To get our estimates we adopt tTinei ~ 130 mb [28]. 
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The inclusive pp ^ G production cross section has been evahiated [29] using CTEQ5L 
parton distribution functions [30]. For Mq = 500 GeV, and integrated over all values of s, 
the cross section can be conveniently parametrized as 



app^G = 1.17 X 10-^^ f-^ ) cm^ . (8) 



1.0565 

VGeV. 

When the condition s > 16Mq is imposed on the integration, we have estimated that the 
cross section given in Eq. (8) is decreased by a factor of 2. The dependence on Mq can be 
roughly described by the scaling behavior (>-0.007iMG-500 GcV)^ 

The Eg integration interval in Eq. (7) is very narrow, as discussed above. Since we 
are working at the very high energy end of the cosmic ray spectrum, we will examine the 
consistency of Eq. (7) with bounds on the neutrino flux at the highest energy. Inserting 
Eq. (8) evaluated at (Ep) = 7.5 x 10^'^ GcV and appropriately reduced to allow for the Smin 
cut, as well as Eq. (4), into Eq. (7), we obtain (with Ng = 2) for the high energy neutrino 
flux accompanying G-production 

ME,) AE, ^ 1200 (N,) J Jg{Eg) dEG (9) 
{N,) 1.7 X 10"^^ cm-^ s~^ sr-^ (10) 

for an event rate of 1 yr~^. At the end of the spectrum, we may approximately consider 
only the most energetic of the secondary tt^, carrying about 8% of the primary energy [31]. 
This entails a neutrino multiphcity {Ni) — 3/event, each carrying 1/4 of the pion energy. 
Thus AE^ ~ {EJ) ~ 0.02 {Ep) = 1.5 x 10^^ GeV. The expected neutrino flux (all flavors) 
at this energy is then 

J^(1.5 X 10^^ GeV) f« 4.4 x 10"^° GeV^^cm'^ s"^ sr"^ . (11) 

An upper bound for the neutrino flux in this energy range has been obtained through the 
absence of radio signals originating in the Moon's rim (GLUE [32]) or in the Greenland 
ice sheet (FORTE satellite [33]). The flux associated with 1 G event/yr at PAG, given in 
Eq. (11), is comfortably below these limits. The RICE experiment [34] (detecting electron 
neutrino-induced radio Cerenkov radiation in the polar ice cap) has reported upper bounds 
on the Ve + i^e flux for energies up to 10^^ GeV. A slight extrapolation of their result to 
1.5 X 10^^ GeV gives a (3-flavor) upper bound which is a factor of 2.5 lower than the flux 
in Eq. (11). The flux in Eq. (11) is also a factor of 3 lower than the model-independent 
bound found in [35] from the absence of horizontal air showers. Finally, extrapolation of 
the ultra-high energy neutrino intensity given in Eq. (11) down to lower energies, assuming 
Ju{Ey) cx E'"^ (more on this below), leads to a neutrino flux which is in agreement with all 
upper limits on J„{Ei,) reported by the AMANDA GoUaboration [36]. 

There is an additional bound (known as the cascade limit) coming from EGRET's observa- 
tion of GeV gamma rays [37]. The relation to neutrinos is as follows: isotopically symmetric 
triplets of n^, , and tt" produced at high energy sources yield 3 i/'s per charged pion and 
2 7's per 7r°, with energies E,^ — E^^/A and Ej — Et^/2, respectively. {Et^ is the pion energy.) 
The integrated neutrino and gamma ray energies then satisfy [38] 



/ E, J,{E,) dE,^- E^ J^{E^) dE^ . (12) 
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Normalization to EGRET data [37], 



oo 

JJE') dE' = 1.45 X 10"^ cm-2 s"^ sr"^ , (13) 

0.1 GeV 

with a spectrum Jf{Ej) — CjE~^, implies C-y = 1.45 x 10~^GeV cm~^ s~^ sr~^ Using 
Eq. (12) we obtain = 3/2C^ for the normalization constant of the accompanying neutrino 
spectrum, again on the assumption of a spectrum J^, oc E^'^. At face value this gives a flux 
at 1.5 X 10^^ GeV which is nearly a factor of 5 smaller than the flux in Eq. (11). However, 
this represents a very large extrapolation: even logarithmic corrections to the E~^ spectrum 
could result in sizeable deviations at very high energies. To see how such corrections might 
arise, we note that for a primary high energy proton with energy Ep, the resulting pion 
spectrum is expected to obey a modified Feynman scaling in the central rapidity region, 
dN„ldETr\Ej, ~ C{Ep) / Ett, where C may be growing as some power of In^Ep) [39]. For 
given Et, < 0.08 -Ep,max, we may convolve with a proton spectrum typical of Fermi engines, 
dNp/dEp oc 1/Ep, to obtain the pion spectrum [40]: 



dN^ fEp,n.s^ dN^ 
dE, 



dE^^ jET^/om ^ dEj^ 



dNp C{E^ (14) 



where C{Et^) is generically a function which grows as a power of \n{ET^), falling to zero at 
the cutoff Et^ = 0.08 ii^p,max- Hence we suspend judgement with respect to the constraint 
imposed by the cascade bound. 



IV. AsuSY WRITTEN IN THE SKY? 

We have described conditions under which G-hadrons could be produced in astrophys- 
ical sources in sufficient abundance to allow detection in an air shower array. Moreover, 
the complete longitudinal and shower-front curvature profiles of G-hadron cascades would 
uniquely differentiate them from proton and nucleus backgrounds. In this last part of the 
paper we examine what we can learn about SUSY if such events are actually observed. 

We have seen that limits on heavy isotope abundance place an upper bound on Asusy- 
Detection of G-hadrons, presumed to originate at cosmological distances, will place a lower 
bound on the proper hfetime Tq of the G. With the use of Eq. (1) this can be translated into 
a lower bound on the scale of SUSY breaking. The argument can be specified as follows. 
In the presence of decay, the integration of the G-flux over source distances r out to the 
horizon 7^ pa 3 Gpc is modified by inclusion of a damping factor 

/ = /^e-"/[^^(^°)ldr (15) 
Jo 

^ (l _ e-^/\criEa)]^^ ^ (16) 

where the Lorentz dilated G lifetime is given in terms of the proper lifetime t{E) = 
TqEq/Mq, with Eq 5 X 10^ GeV. If a few G events are seen during the hfetime of 
PAO, and the observation of > 10^^ GeV neutrinos merits the expectation of 1 G event/yr, 
then / cannot be too small, say / ^ 0.1. This implies that ct{Eg)/TZ ^ 0.1, and places a 
bound 

To > 100 yr (17) 
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on the proper lifetime. Prom Eq.(l), one would then obtain a lower limit on the SUSY 
breaking scale, 

AsusY > 2.6 X 10^^ GeV . (18) 

In conjunction with the upper limit from isotope abundance, Asusy < 10^'^ GeV, we can then 
adduce the remarkable result that the observation of a few G events during the operating 
life of PAO can fix the scale of high energy SUSY breaking. 

V. SUMMARY 

In this work we have examined under which conditions mctastablc G-hadrons can be 
synthesized in powerful astrophysical environments and can be detected on Earth. We 
show, that if cosmic ray sources are able to accelerate protons somewhat above 5x10^^ GeV, 
about 0.5-1 G-hadron induced shower per year could be detected at PAO. Additionally, G- 
cascades provide a signal which is easily differentiated from background. These numbers are 
for Mg = 500 GeV: the event rate could be doubled for Mo = 400 GeV. 

Should some of these very distinctive showers be observed, a lower bound on Asusy can 
be deduced in a manner unavailable at colliders. The combination of such a lower bound 
with the upper bound imposed by the scarcity of heavy isotopes fixes the scale of SUSY 
breaking to a relatively narrow window, 2.6 x 10^^ GcV < Asusy < 10^^ GeV. This result, 
if true, combined with the expected low mass of the higgs would provide strong support for 
a finely-tuned universe. 

Note added: After this paper was completed, a work with related interesting ideas for 
detection of long-lived gluinos at colliders and at IceCubc appeared [41]. For observation at 
IceCube, the gluino production takes place in cosmic ray collisions in the Earth's atmosphere. 
The steeply falling cosmic ray luminosity above 10^ GeV allows sufficient G production for 
a measurable signal at IceCube only for Mq ^150 GcV, which is very close to the present 
bounds from Tevatron [42]. In examining production at extraterrestrial sources, the present 
work can substantially extend the possible range of Mq able to be probed by cosmic ray 
experiments. 
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